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REVIEW  OF  HYDROLOGY  CONCEPTS 


The  science  of  hydrology  has  undergone  enormous  changes  in  the  last 
few  years.    Many  techniques  now  considered  common -pi ace  had  scarcely  been 
conceived  of  ten  years  ago.    Even  the  teminology  has  changed  radically  in 
the  past  decade.    Hydrologists  now  talk  conversationally  of  deterministic , 
parametric,  or  stochastic  hydrology,  of  Markov  chains,  and  Monte  Carlo 
methods.    When  this  research  project  was  undertaken  in  1963  it  would  have 
been  difficult  to  even  find  reference  to  such  terms  in  the  most  advanced 
hydrologic  literature .    In  this  section  an  overlook  at  modern  prediction 
hydrology  will  be  presented  with  an  attempt  to  show  how  the  effects  of  this 
project  relate  to  present-day  hydrologic  practices. 

The  American  Society  of  Civil  Engineers  Committee  on  Surface  Water 
Hydrology  (1963)  described  "parametric  hydrology"  as  the  development 
of  relationships  among  physical  parameters  involved  in  hdyrologic  events  and 
the  use  of  these  relationships  to  generate,  or  synthesize,  non-recorded 
hydrologic  sequences.    The  same  ASCE  committee  described  "stochastic  hydro- 
logy" as  the  use  of  statistical  characteristics  of  hydrologic  variables 
to  solve  hydrologic  problems.    As  related  to  peak  discharge  investigations, 
parametric  hdyrology  methods  all  solve  attempts  to  somehow  describe  a 
"system  mathematically.    The  system  may  include  characteristics  of  a  water- 
shed, meteorological  phenomena  of  a  storm,  etc.    The  system  is  often  pictured 
as  a  "black  box"  into  which  Ideally  one  feeds  a  rainfall  event  (for  example), 
turns  the  crank  and  pulls  out  a  flood  hydrograph.    Stochastic  hydrology 
methods  Is  peak  discharge  studies  are  those  which  pay  no  attention  to 
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precipitation,  watershed  characteristics,  but  rather  consider  only  the 
statistical  nature  of  recorded  runoff  events. 

In  addition  to  prediction  methods  which  are  based  upon  purely  parametric 
hydrology  or  purely  stochastic  hydology,  there  are  a  number  of  techniques 
which  incorporate  some  features  of  each.    An  example  would  be  a  procedure 
in  which  a  number  of  gaged  watersheds  are  studied  by  stochastic  hydrology, 
and  the  results  of  this  study  extrapolated  to  nearby  ungaged  watersheds  on 
the  basis  of  a  regression  equation  relating  watershed  parameters. 

■ 

A  comment  regarding  "deterministic  hydrology"  might  be  useful  at  this 
point.    Deterministic  hydrology  implies  that  a  unique  cause-effect  relation- 
ship exists,  with  no  element  of  chance  involved.    For  a  rainfall -runoff  event 
that  occurred  last  year,  for  example,  it  would  be  theoretically  possible  to 
write  a  singular  mathematical  function  relating  the  streamflow  to  the  causative 
precipitation,  infiltration,  evaporation,  etc.    Hence,  this  is  an  example  of 
deterministic  hydrology.     (The  fact  that  the  mathematical  function  may  never 
be  written  because  of  lack  of  data  does  not  alter  the  idea  that  the  solution 
is  nevertheless  deterministic ) .    With  regard  to  a  rainfall -runoff  event  that 
will  occur  next  year,  it  would  be  impossible  to  define  an  unique  hydrograph, 
perfect  in  every  detail,  for  a  given  rainstorm  (even  if  it  were  possible  to 
know  in  advance  the  precise  characteristics  of  the  rainstorm).    Any  aspect 
of  prediction  hydrology,  therefore,  is  stochastic,  because  there  is  always 
at  least  some  element  of  probability,  chance,  or  statistical  variation 
attached . 
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All  methods  of  peak  discharge  investigation,  whether  parametric, 
stochastic,  or  mixed,  share  two  inherent  limitations,  which  are  not  always 
expressed,  but  which  are  always  implied.    All  depend  upon  historical  records 
of  the  values  of  some  measured  parameters,  and  all  depend  upon  the  assumption 
of  time  invariance  or  "stationarity"  of  the  hydrologic  systems.    The  first 
limitation  implies  that  errors  in  the  basic  measured  data  will  cause  errors 
(perhaps  larger)  in  the  predictions;  the  second  limitation  implies  that  if 
long-term,  natural  or  artificial,  cyclic  changes  or  trends  (the  earth  is  warming 
up,  for  instance)  exist,  errors  will  be  introduced  into  the  studies,  and  no 
prediction  of  the  magnitude  of  these  errors  can  be  made. 

PARAMETRIC  METHODS 

Several  differnt  Parametric  Methods  have  been  proposed  and  these  will 
be  described  herein.     Included  are  Correlation  Analysis,  Combination  System 
Synthesis  and  Analysis,  General  System  Synthesis,  and  General  System  Analysis. 

Although  these  methods  have  certain  basic  differences,  they  all  conceive 
of  the  rainfall -runoff  phenomena  as  a  three-phase  problem.    As  shown  In 
Figure  A-l,  input  (precipitation),  is  acted  upon  by  a  system,  to  produce  an 
output  (runoff) . 

Given  two  of  the  three  phases  shown,  the  third  can  be  obtained  (at  least 
a  solution,  but  not  necessarily  an  unique  solution).    Amorocho  and  Hart 
pointed  out  (196-4)  that  in  the  usual  investigation,  recorded 
historical  Input  and  output  data  are  used  (with  or  without  watershed  para- 
meters), to  construct  a  mathematical  model  of  a  system  which  (hopefully) 
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accurately,  describes  the  real-world  system.    This  mathematical  model  is 
then  used  in  conjunction  with  hypothesized  (synthetic)  input  or  output  to 
generate  predicted  values  for  the  unknown  phase*    Operation  of  the  model 
after  it  has  been  constructed  (which  of  course  is  the  purpose  or  motive 
for  the  investigation  in  the  first  place)  is  abbreviated  by  Ameroeh©  and  Hart 
as  BQF  (for  reconstruction  of  past  sequences  or  prediction  of  future  behavior). 

an  inp1 

r©  c*^i.Y*  jr©s^      ^i*©  c^xiq  n  g  i  ©  s         ^^^^^^  $ 
Correlation  Analysis 

in  this  proeeture  various  oojMnations  of  variables  are  tested  to  explore 
tfe§  significance  of  tfeeir  tffeets  in  the  hyirol©§i§  system*   For  a  §£yb& 

combination  of  variables  a  regression  analysis  is  performed  to  determine  the 

best^fit  equation  relating  %km%  Satfe  mok  eausMon  (or  motel)  yieiis  s  sot 
of  error  e§tiia.t©§  §n.i  measufej  §f  eemlsliea  i  m  ittepiinii  W  sutlysis 
©f  var4sn§§*   ffiii  §©mMn§M§!i  of  variables  wkieh  yieiis  i  relationship  most 
elosely  spiPrtiating  tie  r©c©ri§ i  output  fwttfeloi  in  tens  of  %)m  reetHsi 

iBfttt  f^lgtion  tm  determine^  by  the  error  estimates  end  measures  of  oorrola-? 

ti§n)  4§  £§§!&§£  I©  tee  %lm  tees*  ftiA4ftt&4§R.j  sni  Its  r@gieg§i§ii  ©tuaMoi  is 

adopted  as  the  prediction  equation. 


\ 


The  often-used  and  much-abused  "rational  formula" 

Q  =  C  i  A  (A_!) 

might  be  visualized  as  a  correlation  equation  although  it  is  doubtful  whether 
it  was  originally  formulated  as  the  result  of  regression  analysis.    It  relates 
peak  discharge  Q  to  an  input  term  i  (a  statistical  value  of  average  rainfall 
intensity  over  a  time  period  corresponding  to  the  estimated  time  of  concentration 
of  the  drainage  area)  and  a  watershed  parameter  A  (watershed  area) .  The 
runoff  coefficient  C  could  be  thought  of  as  a  regression  coefficient  or 
retardance  coefficient  which  has  been  determined  empirically  for  various 
ground  cover  conditions.    Many  equations  more  sophisticated  than  the  rational 
formula  have  been  proposed,  involving  linear  or  nonlinear  relationships 
among  few  or  many  variables . 

Correlation  analysis  methods  are  basically  empirical.    Their  principal 
disadvantages  are  that  almost  equally  good  prediction  equations  often  result 
from  entirely  different  sets  of  parameters,  and  the  investigator's  judgment 
as  to  the  physically  relevant  parameters  largely  controls  the  success  of  the 
operation.    Also,  use  of  parameters  which  are  interdependent  may  lead  to 
erroneous  results.    Sharp,  et  al,  (i960)  reported  on  multiple  regression 
studies  of  two  watersheds,  one  in  Kansas  and  one  in  Oklahoma.    They  used 
various  combinations  of  up  to  eight  independent  parameters  and  one  dependent 
parameter  (discharge  peak).    They  found  that  the  quality  of  the  predicting 
equations  became  progressively  poorer  for  values  far  removed  from  the  average. 
Although  they  had  high  correlation  coefficients  this  seemed  to  indicate  only 
that  reasonably  good  individual  estimates  could  be  made  near  the  means,  but 
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as  estimates  were  made  for  values  farther  and  farther  removed  from  the  means , 
probable  errors  increased  enormously.    In  many  cases  they  found  inverse 
correlations  with  some  of  the  variables  when  physical  considerations  plainly 
seemed  to  require  direct  correlations,  and  vice  versa.     They  concluded  that 
their  indifferent  success  could  be  traced  to  four  assumptions  inherent  in 
application  of  the  methods: 

1.  There  are  errors  only  in  the  dependent  variable  (discharge  peak). 
No  errors  occur  in  the  independent  variables. 

2.  The  variance  of  the  dependent  variable  does  not  depend  on  the 
values  of  the  independent  variables. 

3.  The  observed  values  of  runoff  are  uncorrelated  random  variables. 
h.      The  population  of  the  dependent  variable  is  normally  distributed 

about  the  regression  line  for  any  fixed  level  of  the  independent  variables 
(precipitation ,  for  instance). 

Sharp,  et.  al. ,  felt  that  with  hydro logic  data  all  these  assumptions 
are  suspect  and  that  the  fourth  assumption  in  particular  is  usually  violated 
because  the  discharge  peaks  are  badly  skewed  about  the  regression  line  and 
do  not  approach  a  normal  distribution.     For  these  reasons,  and  because  there 
seems  to  be  no  reliable,,  nonsub jective  way  to  eliminate  personal  factors 
associated  with  the  selection  of  relevant  parameters  to  begin  with,  they  con- 
cluded that  very  little  reliance  could  be  placed  on  the  final  equations. 
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Combined  System  Analysis  and  Synthesis 


To  understand  this  method  and  those  that  follow,  it  is  first  necessary 
to  define  and  discuss  the  terms  "system  analysis"  and  "system  synthesis." 
In  system  analysis  a  mathematical  process  involving  the  use  of  measured 
input  and  output  only  (without  any  attempt  to  describe  explicitly  the 
internal  mechanisms  of  the  system)  is  used  to  establish  the  relationship 
between  input  and  output.     This  relationship  takes  the  form  of  a  unique 
function  which  is  made  to  operate  on  the  input  to  produce  the  output.  The 
function  does  not  have  to  have  a  physical  meaning  or  possess  parameters 
fulfilling  conditions  of  dimensional  consistency.    The  analysis  may  be  either 
linear  or  nonlinear  depending  upon  the  form  of  the  function  selected  to 
establish  the  input-output  relationship. 

In  system  synthesis  the  operation  of  the  system  is  described  by  a  linkage 
or  combination  of  components  whose  presence  is  presumed  to  exist  in  the  system 
and  whose  functions  are  known  and  predictable.     The  linkage  of  components  is 
made  in  such  a  manner  that  whenever  a  specified  input  is  applied  the  correct 
output  is  produced;  but  the  process  of  synthesis  usually  does  not  yield  an 
unique  model  of  the  system. 

Pure  analysis  or  pure  synthesis  can  be  performed  on  a  system  independently, 
but  frequently  a  combination  of  both  is  employed.     The  classical  unit 
hydrograph  procedure  is  an  example  of  a  combination  of  analysis    and  synthesis . 
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Unit  Hydrograph 

The  unit  hydrograph  procedure  is  a  three-phase  operation  as  has  been 
described  and  diagrammed  earlier  (page  A- 3  ).     In  this  case  the  gross  input 
in  first  modified  by  infiltration  over  the  surface  of  the  watershed  to  yield 
a  "rainfall  excess"  function.    The  surface  elements  of  the  watershed  then 
operate  in  the  rainfall  excess  function  to  generate  a  surface  runoff  function, 
which  is  modified  in  turn  by  the  base  flow  to  produce  a  flood  hydrograph. 
The  first  and  third  phases  of  the  procedure  could  be  classed  as  system 
synthesis,  while  the  second  (inversion  of  rainfall  excess  into  surface  runoff 
by  means  of  a  unit  hydrograph)  is  an  example  of  system  analysis. 

The  first  phase,  modification  of  gross  input  into  rainfall  excess 
function,  is  accomplished  by  one  of  several  techniques  each  of  which  is 
empirical,  subjective,  and  inexact.     Infiltration  tests  might  be  run  on 
different  soils  on  the  watershed  to  determine  infiltration  capacities  (maxi- 
mum rates  at  which  infiltration  can  occur);  soil  moisture  conditions  prior 
to  the  storm  might  be  estimated  (by  actual  soil  moisture  measurements  if 
possible,  or  by  use  of  an  'antecedent  precipitation  index*);  or  judgment 
might  be  used  to  determine  that  some  fraction  of  gross  rainfall  should  be 
expected  to  infiltrate  into  the  soil  or  otherwise  be  lost  through  surface 
detention,  interception,  etc.    The  resulting  rainfall  excess  function 
is  certainly  not  determined  uniquely,  but  it  does  result  from  attempts  to 
describe  the  operation  of  this  part  of  the  system  by  means  of  components 
which  are  undeniably  present. 

The  third  phase,  that  of  modifying  the  surface  runoff  function  into  flood 
hydrograph,  is  again  an  inexact  procedure.    There  is  no  doubt  as  to  the 
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existence  of  base  flow,  nor  the  fact  that  it  modifies  storm  runoff.  But 
subjective  techniques  must  be  employed  to  make  the  actual  transformation 
and  the  resulting  solution  is  not  unique. 

The  second  phase,  which  of  course  is  the  heart  of  the  procedure, 
employs  the  unit  hydrograph  concept  advanced  by  Sherman  (1932).  Sherman 
defined  the  unit  hydrograph  as  that  hydrograph  resulting  from  one  inch  of 
rainfall  excess  from  a  storm  of  specified  duration.    Since  losses  and  accretions 
are  accounted  for  in  phases  1  and  3,  the  unit  hydrograph  represents  exactly 
the  same  volume  of  water  as  did  the  rainfall  excess  hyetograph,  namely  one 
inch  of  water  from  the  entire  watershed  area.    It  therefore  represents  a 
direct  inversion  of  water  from  rainfall  to  streamflow. 

A  unit  hydrograph  is  developed  for  a  watershed  by  analysis  of  a  suitable 
storm  hydrograph,  after  first  subtracting  out  the  base  flow.    Knowing  the 
volume  of  runoff  in  the  surface  runoff  hydrograph,  each  ordinate  is  scaled 
upward  or  downward  a  proportionate  amount  to  produce  a  hydrograph  having  the 
required  one  inch  of  volume.    This  is  said  to  be  the  unit  hydrograph  for 
storms  of  similar  duration.    After  it  has  once  been  developed,  the  unit 
hydrograph  may  be  used  for  ROP  purposes  and  is  strictly  a  linear  analysis 
tool.     (A  storm  with  two  inches  of  rainfall  excess  will  result  in  a  surface 
runoff  hydrograph  with  a  peak  discharge  exactly  double  that  from  a  storm 
with  one  inch  of  rainfall  excess).    It  gives  unique  results  which  are  not 
dependent  on  watershed  factors,  etc. 
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In  practice  several  storm  hydro graphs  (from  storms  of  similar  duration) 
may  be  analyzed  to  produce  an  average  unit  hydrograph,  and  of  course  many 
unit  hydrographs  are  needed,  one  for  each  storm  duration.     (Linsley ,  Kohler, 
and  Paulhus ,  however,  suggested  (1958)  a  method  whereby  available  unit  graphs 
for  a  watershed  might  be  utilized  in  producting  a  unit  graph  for  some  other 
storm  duration). 

On  small  watersheds,  and  in  arid  regions  generally,  data  are  not  often 
available  from  which  to  construct  unit  hydrographs.    Snyder  (1938)  pro- 
posed a  "synthetic  unit  hydrograph"  method  for  use  in  such  situations.  He 
defined  the  unit  hydrograph  peak  in  terms  of  drainage  area  and  time  lag 
(the  time  interval  between  centroid  of  rainfall  excess  and  the  hydrograph 
peak).    Another  empirical  relation  is  used  to  define  the  shape  of  the  hydro- 
graph  ;  and  the  unit  hydrograph  thus  synthesized  may  be  used  in  the  same 
fashion  as  ordinary  unit  hydrographs.    Unit  hydrograph  synthesis  can  only 
be  applied  to  watersheds  where  major  floods  result  from  single  storms  over 
the  entire  basin,  and  this  restriction  limits  the  procedure  to  use  on  very 
small  basins  only.    Hickock,  et.al. ,  (1959)  suggested  that  in  the 
Southwest  synthetic  hydrographs  could  not  be  used  on  watersheds  larger 
than  1000  acres.    Renard  and  Keppel  ( 1966 )  reported  that  they  had  divided 
larger  watersheds  into  sub are as ,  developed  synthetic  unit  hydrographs  for 
each,  and  had  used  routing  techniques  to  obtain  final  hydrographs  at  water- 
shed mouth.    These  methods  are  extremely  time  consuming  and  laborious. 
The  U.  S.  Soil  Conservation  Service  uses  essentially  the  same  procedure  for 
ungaged  watersheds,  but  utilizes  digital  computer  programs  to  expedite 
solution.     The  SCS  technique  was  described  in  detail  by  Ferris  and  Williams 
(1968)  and  has  been  discussed  earlier  in  this  report   (Chapter  VI ). 
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The  unit  hydrograph  procedure  has  been  widely  acclaimed  as  a  valuable 
tool  in  flood  predicting.    Recently,  however,  there  have  been  several  disturb- 
ing reports  which  imply  that  the  procedure  is  less  valid  than  had  been  though. 
Ancorocho  and  Hart,  for  instance,  cite  evidence  indicating  that  this 
method  tends  to  underestimate  flood  episodes  which  are  larger  than  those  used 
for  its  analysis  and  to  overestimate  small  floods .    They  attribute  this 
apparent  defect  to  the  fact  that  the  watershed  is  basically  a  nonlinear 
system,  to  which  is  forced  a  linear  function.    They  concede  that  there  are 
many  instances  where  linear  approximations  are  valuable,  and  that  linear 
analysis  is  much  simpler  and  easier  to  use  than  non-linear  analysis,  but 
they  caution  that  grossly  erroneous  predictions  may  accrue,  without  any 
foreknowledge  of  the  magnitude  or  sign  of  the  error.    These  authors  conclude 
by  pointing  out  that  the  unit  hydrograph  method  treats  the  input  and  output 
as  'lumped'  parameters,  meaning  that  the  areal  distribution  and  variability 
of  the  input  (rainfall)  are  not  considered  to  have  any  effect  on  the  output. 
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General  System  Synthesis 


With  the  advent  of  the  electronic  computer  has  come  the  capability  for 
construction  of  synthetic  hydrologi c  watershed  models,    The  much-publicized 
Stanford  Watershed  Model  (Crawford  &  Linsley ,  19 66)  is  the  best  known  example 
of  this  method.     The  Stanford  Model  (and  others  like  it)  rely  on  a  statement 
of  continuity,  which  can  be  expressed  as 

I  =  Q  +  AS  (A-2) 

where  I  is  total  inflow,  Q  is  total  outflow,  and  AS  is  change  in  internal 
storage;  all  referred  to  a  time  interval  t^  to  t^.    This  equation  is 
differentiated  to  yield 

i(t)  =  q(t)  +  (—-)  (A-3) 

where  i(t)  is  rate  of  total  inflow,  q(t )  is  rate  of  total  outflow,  and 
ds/dt  is  rate  of  change  in  internal  storage. 

The  process  of  general  synthesis  typically  is  begun  by  postulating 
a  simulation  model  (usually  complex)  whose  structure  is  based  on  an  assumed 
knowledge  of  the  phenomena  involved  in  the  hydrologic  cycle.    The  Stanford 
Model  for  instance  includes  relationships  involving  precipitation,  evapo- 
transpiration ,  interception ,  infiltration ,  groundwater  storage  and  movement , 
overland  flow,  interflow  and  channel  delay  (among  others).    Explicit  mathe- 
matical functions  are  developed  to  define  operations  effected  on  various 
elements  of  the  input  and  storage ,  after  which  the  recorded  input  is 
processed  through  the  model ,  usually  with  the  aid  of  a  digital  computer 
(occasionally  an  analog  computer)  and  the  resulting  output  compared  with 
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recorded  output  of  the  natural  system.     If  there  is  unsatisfactory  agree- 
ment, various  components  of  the  model  are  adjusted  or  modified,  and  the 
process  repeated  until  satisfactory  agreement  between  computed  output  and 
recorded  output  results.    At  that  point  the  model  is  assumed  to  be  properly 
adjusted  and  it  may  then  be  used  for  ROP  purposes. 

Most  of  the  synthetic  models  developed  thus  far,  including  the  Stanford 
model,  have  proceeded  from  a  very  detailed  model  of  a  relatively  small  area. 
As  computer  capabilities  and  modelling  experience  have  progressed,  some  of 
these  models  have  been  expanded  to  encompass  larger  and  larger  areas.  An 
investigation  which  has  been  recently  undertaken  by  the  departments  of  Civil 
Engineering  and  Engineering  Mechanics  and  Industrial /Management  Engineering 
at  Montana  State  University  is  aimed  at  development  of  a  state-wide  model  for 
Montana.     This  model  is  based  on  a  "systems"  approach  of  beginning  with  a 
state  boundary  model,  progressively  subdividing  into  smaller  and  smaller  units, 
with  corresponding  increases  in  completely  defining  system  components  in 
greater  detail  and  with  greater  accuracy. 

Both  linear  and  nonlinear  models  have  been  produced.    As  was  mentioned 
earlier,  linear  relationships  are  far  easier  to  develop  than  are  non-linear 
relationships.    But  ultimately  the  basic  non-linearities  in  the  basin  must  be 
reckoned  with  and  recognized  in  the  model. 

Reports  in  the  literature  indicate  that  results  from  some  of  the  water- 
shed models  are  encouraging;  but  there  are  differences  of  opinion  as  to 
how  reliable  they  are  for  reconstructing  or  predicting  detailed  runoff  events, 
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such  as  peak  flow  hydrographs .    Amorocho  and  Hart  suggest  four  probable  causes 
for  lack  of  reliability:     errors  in  recorded  data  from  which  the  model  is 
constructed,  effects  of  areal  distributions  of  parameters,  imperfections 
of  the  model  structure,  and  nonuniqueness  of  the  synthesis  process. 

The  first  two  of  these  causes  are  common  to  other  techniques  as  well 
as  to  that  of  system  synthesis,  but  these  authors  fear  the  magnitude  of 
errors  may  be  much  larger  when  using  synthesis.    Subdivision  of  the  hydro- 
logic  unit  into  subareas ,  each  with  its  own  parameters,  is  a  way  of  reducing 
the  problem  of  are ally  distributed  parameters.    This  technique  is  used  on 
the  Stanford  model,  and  is  being  built  into  the  Montana  Water  Planning  model. 
Model  structure  imperfections  are  inevitable  because  of  the  necessity  to 
simplify  the  model  from  the  natural  system.    In  the  first  place  a  natural 
system  is  simply  too  complex  to  describe  in  all  its  details;  in  the  second 
place,  the  kinds  and  amounts  of  data  required  would  become  nearly  impossible 
to  collect.    Since  the  synthesis  process  is  not  unique,  it  follows  that  any 
synthetic  model  must  of  necessity  be  non-exact.    Even  though  a  model  faith- 
fully reproduces  historic  records,  there  must  always  be  uncertainty  as.  to 
its  performance  when  out-of- range . 
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General  System  Analysis 

Most  system  analysis  attempts,  at  least  up  until  a  few  years  ago,  have 
been  based  on  linear  techniques.    The  unit  hydrograph  method,  for  example, 
utilizes  a  linear  analysis  (with  synthesis  involved  at  both  input  and  out- 
put ends).    As  was  mentioned  earlier,  the  watershed  is  basically  a  nonlinear 
system,  and  the  unit  hydrograph  procedure  of  forcing  all  the  nonlinear! ty 
into  the  input  and  output  synthesis  is  artificial.     Linear  techniques  are 
much  preferred  by  investigators  because  of  their  simplicity.    In  fact,  the 
analysis  of  nonlinear  systems  had  been  barely  considered  prior  to  I960. 
Within  the  past  few  years,  however,  considerable  progress  has  been  made  in 
the  adapting  of  nonlinear  analysis  to  the  study  of  hydrologic  systems . 

An  advantage  to  the  use  of  general  system  analysis  is  the  fact  that  it 
is  the  fact  that  it  permits  relationships  to  be  established  between  partial 
inputs  and  partial  outputs  without  satisfying  physical  continuity  conditions 
for  the  entire  system.    In  theory  at  least  it  permits  direct  operation  with 
recorded  precipitation  and  recorded  river  runoff  without  considering  internal 
parameters  such  as  evapotranspiration  and  groundwater  flow. 

The  process  consists  essentially  of  analyzing  the  recorded  input  and 
recorded  output  to  find  an  equivalent  system  function  or  convolution.  In 
the  linear  case  the  inversion  is  the  unit  hydrograph,  or  its  equivalent; 
in  the  nonlinear  case  the  inversion  is  an  integral  function. 

Current  procedures  utilizing  nonlinear  analysis  usually  involve  approxi- 
mations because  general  methods  of  direct  nonlinear  inversion  are  not  yet 
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available.    Araorocho  and  Orlob  (1961)  reported  a  nonlinear  system  analysis 
which  was  based  on  the  assumption  that  a  hydrologic  system  is  equivalent  to 
an  analytic  system  of  progressively  higher  ordered  terms,  whose  output  when 
added  equal  the  output  of  the  natural  system; 


y(t)  =  hx  (t^  x  (t  -  rx)  d  z~1 


♦  r 


OO  J 


\  (r^rj  x  (t  -rx)  x  (t  -  r2)dr1  ar2  (A"M 


+  r  /y  . . . » +  


where  h^,  h^,  etc.  are  system  functions  of  progressively  higher  order;  x  is 
input;  y  is  output.    All  terms  are  functions  of  time.    The  task  is  to 
evaluate  each  of  the  nonlinear  kernels  of  the  functions  of  this  series.  When 
complex  inputs  are  involved,  approximation  is  necessary  to  evaluate  the  series 


In  its  present  stage  of  development,  nonlinear  analysis  has  been 
applied  only  to  systems  which  are  time  invariant,  and  to  systems  with  lumped 
parameters.    The  method  offers  promise,  however,  of  utility  in  distributed 
parameter  situations  with  time  variance. 


-  A- 17  - 


STOCHASTIC  METHODS 


The  parametric  methods  described  above  for  determining  streamf low 
discharges  are  not  often  satisfactory  for  establishing  return  frequency 
directly.     Usually  these  methods  require  starting  from  a  hypothesized  rainfall 
event,  generating  a  hydrograph ,  and  finally  estimating  the  return  frequency 
of  the  causal  storm. 

Stochastic  hydrologic  methods  are  used  to  establish  the  level  of 
probability  with  which  certain  streamf low  magnitudes  or  sequences  of  magnitudes 
are  likely  to  occur  in  the  future.    They  work  strictly  from  previously 
recorded  streamf low  data,  and  are  in  no  way  dependent  upon  watershed  or 
storm  event  parameters.    The  stochastic  problem  is  simply  this:     given  a 
historic  sequence  of  events,  with  historic  data  available,  what  inferences 
can  be  deduced  from  their  statistical  distribution  to  assess  the  probabilities 
of  future  sequences . 

Mass  Diagrams 

Until  recent  years  analysis  of  reservoir  capacity  requirements,  possible 
yield,  and  other  similar  problems  usually  started  with  the  assumption  that 
the  past  recorded  record  fully  describes  all  possible  combinations  of  events, 
and  that  future  operations  can  be  planned  by  utilizing  such  a  chronology 
without  transformation.     In  spite  of  the  weaknesses  of  such  an  assumption 
this  procedure  is  still  widely  used  for  these  purposes.     Construction  of  a 
mass  diagram,  which  shows  cumulative  flow  volume  as  a  function  of  time, 
permits  graphical  solution  of  the  problems  being  studied. 


-  A- 18  - 


Statistical  Methods  for  Estimating  Frequencies 

The  traditional  approach  for  determining  peak  discharge  frequencies  by 
stochastic  analysis  has  been  to  conduct  a  frequency  analysis  in  the  light 
of  a  pre -conceived  distribution  model.     At  least  a  dozen  distribution  models 
have  been  proposed,  and  one  of  the  investigator's  biggest  problems  is  deciding 
which  of  these  is  the  correct  one  to  use.    Although  there  is  no  conclusive 
evidence  at  present  to  indicate  that  any  one  model  is  proper  for  all  hydrologic 
regions,  the  federal  agencies  dealing  in  water  resource  development  have  all 
been  directed  by  the  Hydrology  Committee,  Water  Resources  Council  (1967)  to 
use  a  single  method,  the  Log-Pearson  Type  III.    This  model  and  others  will 
be  discussed  below,  following  a  discussion  of  central  tendency  and  divergence. 

Mean  and  Standard  Deviation:     Several  terms  are  used  to  describe  an 
average  or  central  value  of  a  series  of  values  (say  streamflow  peak  discharges). 
In  hydrologic  analyses,  the  measure  customarily  used  is  the  mean ,  which  is 
defined  as 


(A-5) 


where yu  is  population  mean 

is  magnitude  of  an  individual  event 
N  is  the  number  of  events 

Likewise ,  there  are  various  ways  of  describing  the  divergence  or  scatter 
of  events  about  the  average  or  mean.     In  general  no  single  term  gives  a 
complete  description  of  the  divergence,  but  the  most  commonly  used  quantity 
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is  standard  deviation, 


2 


(T     ■  i        1    '  (A-6) 
N  =  1 

where  6^  is  standard  deviation  and  the  other  terms  are  as  defined  before. 


Population  versus  Sample  Mean  and  Standard  Deviation:     In  hydrological 
analyses  "population"  mean  and  standard  deviation  can  not  be  calculated 
precisely  because  the  population  of  hydrologic  events  is  infinite,  and 
values  for  only  a  fraction  of  the  events  are  recorded  or  "sampled."  The 
"sample"  mean  and  "sample"  standard  deviation  are  computed  from  n  recorded 
events  as 


2x. 

X  =  (A-T) 


and 

=T  VA.      -     A.  J  (A_8) 


.2     2T(X.   -  X)2 


n  -  1 

where  X  and  S  are  sample  mean  and  sample  standard  deviation.     It  is  possible, 
of  course,  that  X  and  S  will  have  the  same  numerical  values  asyLx.  and  5"' 
respectively,  but  it  is  more  likely  that  the  values  will  not  be  the  same. 
As  n  becomes  larger,  the  computed  values  of  X  and  S  can  be  expected  to  more 
closely  approximate  .u.  and  (P.     Statistical  procedures  are  available  for 
estimating  the  likely  error  in  X  and  S  for  given  sets  of  data,  and  hence 
setting  upper  and  lower  probable  limits  for    lk  and  (T. 

Probability:     The  likelihood  that  a  certain  event  will  occur  may  be 
expressed  in  terms  of  its  "probability."    Probability  is  a  fraction  lying 
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between  0  and  1;  the  larger  the  fraction,  the  greater  is  the  likelihood  that 
the  event  does  occur.     If  the  probability  is  1,  it  is  mathematically  certain 
that  the  event  will  occur;  if  it  is  zero,  it  is  just  as  certain  that  the 
event  will  not  occur.     Probability  is  defined  by 

p  _   Number  of  successes   , 

Total  number  of  events  (equal  to  number  of  successes  plus  number  of 
failures ) 

If  the  event  being  described  does  occur,  it  is  considered  a  "success";  if  it 
does  not  occur,  it  is  considered  a  "failure".    In  terms  of  streamf low  discharge 
on  a  given  river,  for  example,  where  historic  records  are  being  studied,  the 
probability  that  the  discharge  will  exceed  a  specified  magnitude  could  be 
determined  from 

P  -  Number  Qf  recorded  streamf low  events  greater  than  specified  magnitude 
Total  number  of  recorded  streamf low  events  that  have  occurred. 


Normal  Distribution:     The  equation  of  a  normal  curve  was  published  as 
early  as  1733  by  de  Moivre,  and  developed  independently  by  Laplace,  17^9-1827 , 
and  Gauss,  1777-1855  •    The  normal  curve  was  largely  ignored  until  192k  when 
Karl  Pearson  found  a  copy  of  de  Moivre 's  report  in  a  library.     If  a  series  of 
streamflow  peaks  fits  a  normal  distribution,  the  graph  of  probability  of 
occurrence  versus  magnitude  of  events  produces  a  "bell-shaped"  curve  which 
is  symmetrical  about  the  mean. 

A  typical  normal  curve  (sometimes  referred  to  as  Laplacian  or  Gaussian) 

is  shown  in  Fig.  A-2.  The  curve  theoretically  extends  to  infinity  in  both  direction 
indicating  there  is  no  magnitude  of  Q  for  which  the  probability  of  occurrence 
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is  exactly  zero.    The  area  under  the  curve  between  two  values  of  Q,  say  Q    and  Q 

1  ! 

represents  the  probability  of  the  streamf low  peak  lying  between  those  values 
of  Q.     (This  is  shown  as  the  cross -hatched  area  on  the  figure).    The  total 
area  under  the  curve  is  1,  since  it  represents  all  possible  events. 

For  a  normal  distribution,  the  curve  is  completely  defined  (and  hence 
the  probability  of  any  discharge  occurring  can  be  determined)  if  the  mean, 
t\,  and  standard  deviation,^  are  known.     The  standard  deviation  defines 
a  discharge,  measured  in  each  direction  from  the  mean,  such  that  the  area 
under  the  probability  curve  between  the  discharges      _(J~  and  yiA+T*  is  0.68. 
This  means  that  for  a  normal  distribution  the  probability  is  0.68  that  a 
given  streamflow  peak  will  lie  between  the  values  aa-(T  and  yA+(T. 

The  probability  of  occurrence  of  any  discharge  Q  can  be  computed  from 
Q  =  ^+  Z<T.  (A-10) 
The  multiplier  Z  is  a  frequency  factor  which  can  be  determined  for  any  desired 
probability,  from  Table  A-I .    The  tabulation  simply  describes  the  shape 
of  the  frequency  distribution  function  for  the  case  where  *a  =  0  and  J^=  1.0. 

Because  of  its  simplicity,  the  normal  distribution  is  often  assumed  in 
hydrologic  studies.    But  if  the  normal  distribution  were  correct  it  would 
require  the  possibility  of  negative  discharges,  as  indicated  in  the  sketch 
below.     Sometimes  the  theoretical  probability  of  negative  Q  values  is  small 
enough  so  this  defect  can  be  neglected.     It  is  evident,  however,  that  stream- 
flow  data  never  fit  a  normal  distribution  exactly. 


-  A-22  - 


TABLE  A-I:     Probability  of  a  Random  Value  of    z    =     (  X  -  jla  )/£T  Being  Greater  Than 

the  Values  Tabulated  in  the  Margins  ' 
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.0001 

.0001 

.0001 

.0001 

.0001 

3.9 

.0000 
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Log-Normal  Distribution;     Chow  (195*0  noted  that  many  kinds  of  statistical 
data,  including  streamflow  data,  form  a  skew  frequency  distribution  as  sketched  in 

Fig/  .  A-'S.j'  bu't  that  their _ logarithms  are  nearly: -normally, distributed.  '  If  the 
log  -  normal  distribution  applies,  then  the  discharge  corresponding  to  any 
particular  probability  is  determined  from 

Q  =  y>+  K  <?  (A-ll) 

As  in  the  normal  distribution,  K  is  a  frequency  factor,  but  it  now  depends  on 
the  coefficient  of  skew  or  coefficient  of  variation  as  well  as  probability. 
The  coefficient  of  variation  is 


C    =    -  (A-12) 

V  iJL 


and  coefficient  of  skew  is 


n     £(*3)  -    3  X  2{x2)  +  2  (x) 


C    =  -~   _n  n  (A-13) 

s      n-1         "     —  —  — — — ■ — — — 


Values  of  K  were  developed  by  Chow  from  theoretical  considerations  of  the 

normal  distribution  of  logarithms .     Chow  tabulated  values  of  K  for  selected 

probabilities  and  selected  skew  coefficients.    This  tabulation  is  reproduced 

as  Table  A-II.    The  extreme  right  hand  column  in  this  tables  gives  values 

of  C    which,  in  connection  with  the  coefficient  of  skew  in  the  first  column 
v 

will  produce  plotting  points  of  a  straight  line  on  log-probability  paper. 
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FIGURE  A- 3.     Example  of  a  skew- frequency  distribution. 
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Coeffi-  Proba- 
cient  of  bility 
Skew      at  Mean 


Probability  in  Percentage 
Equal  to  or  Greater  Than  the  Given  Variate 


99       95       80  50 


0 

0.1 
0.2 
0.3 

0.4 

0.5 
0.6 

0.7 
0.8 

0.9 

1.0 

1. 1 

1.2 
1.3 
1.4 

1.5 
1.6 

1.7 
1.8 

1.9 

2.0 
2.1 
2.2 
2.3 
2.4 

2.5 
2.6 

2.7 
2.8 

2.9 

3.0 
3.2 
3.4 
3.6 

3.8 

4.0 

k.5 

5.0 


50.0 
49-3 
48.7 
48.0 

47.3 

46.7 
46.1 

45.5 
44.9 

44.2 

43.7 
43.2 

42.7 
42.2 
4l.7 

4l.3 
40.8 
4o.4 
4o.o 
39.6 

39.2 
38.8 
38.4 
38.1 
37.7 

37.4 

37.1 
36.8 
36.6 
36.3 

36.0 
35.5 
35.1 
34.7 
34.2 

33.9 
33.0 

32.3 


2.33 
2.25 
2.18 
2.11 
2.04 

1.98 
1.91 
1.85 
1.79 
1.74 

1.68 
1.63 
1.58 
1.54 
1.49 

1.45 
l.4i 
1.38 
1.34 
1.31 

1.28 

1.25 
1.22 
1.20 
1.17 

1.15 
1.12 
1.10 
1.08 
1.06 

1.6k 
1.01 
0.98 

0.95 
0.92 

0.90 
0.84 
0.80 


1.65 
1.62 

1.59 
1.56 
1.53 

1.49 
1.46 

1.43 
1.40 

1.37 

1.34 
1.31 
1.29 
1.26 
1.23 

1.21 
1.18 
1.16 
1.1k 
1.12 

1.10 
1.08 
1.06 
1.04 
1.02 

1.00 
O.99 

0.97 
O.96 

0.95 

0.93 
0.90 
0.88 
0.86 
0.84 

0.82 
0.78 
0.74 


0.84 
O.85 
0.85 
0.85 
0.85 

0.86 
0.85 
0.85 
0.84 
0.84 

0.84 
0.83 
0.82 
0.82 
0.81 

0.81 
0.80 

0.79 
O.78 
O.78 

0.77 
0.76 
O.76 

0.75 
0.74 

0.74 

0.73 
0.72 
O.72 
0.71 

0.71 
O.69 
0.68 
0.67 
0.66 

0.65 
0.63 
0.62 


0 

0.02 
o.o4 
0.06 
0.07 

0.09 
0.10 
0.11 
0.13 
0.i4 

0.15 
0.16 
0.17 
0.18 
0.19 

0.20 
0.21 
0.22 
0.22 
0.23 

0.24 
0.24 

0.25 
0.25 
0.26 

0.26 
0.26 
0.27 
0.27 
0.27 

0.28 
0.28 
0.29 
0.29 
0.29 

0.29 
0.30 
0.30 


20          5  1        0.1  0.01 

+         +  +         +  + 

0.84  1.64  2.33  3.09  3.72 

0.84  1.67  2.4o  3.22  3-95 

0.83  1.70  2.47  3.39  4.18 

0.82  1.72  2.55  3.56  4.42 

0.81  1.75  2.62  3.72  4.70 

0.80  1.77  2.70  3.88  4.96 

0.79  1.79  2.77  4.05  5.24 

0.78  1.81  2.84  4.21  5.52 

0.77  1.82  2.90  4.37  5.8i 

O.76  1.84  2.97  4.55  6.11 

0.75  1.85  3.03  4.72  6.40 

0.73  1.86  3.09  4.87  6.71 

0.72  1.87  3.15  5.04  7.02 

0.71  1.88  3-21  5.19  7.31 

0.69  1.88  3.26  5.35  7.62 

0.68  1.89  3.31  5.51  7.92 

0.67  1.89  3.36  5.66  8.26 

0.65  1.89  3.40  5.80  8.58 

0.64  1.89  3.44  5.96  8.88 

0.63  1.89  3-48  6.10  9.20 

0.61  1.89  3.52  6.25  9.51 

0.60  1.89  3-55  6.39  9-79 

0.59  1.89  3.59  6.51  10.12 

0.58  1.88  3.62  6.65  10.43 

0.57  1.88  3.65  6.77  10.72 

0.56  1.88  3.67  6.90  10.95 

0.55  1.87  3.70  7.02  11.25 

0.54  1.87  3.72  7.13  11.55 

0.53  1.86  3.74  7.25  11.80 

0.52  1.86  3.76  7.36  12.10 

0.51  1.85  3.78  7.47  12.36 

0.49  1.84  3.81  7.65  12.85 

0.47  1.83  3.84  7.84  13.36 

0.46  1.81  3.87  8.00  13.83 

0.44  1.80  3.89  8.16  14.23 

0.42  1.78  3.91  8.30  14.70 

0.39  1.75  3.93  8.60  15.62 

0.37  1.71  3-95  8.86  16.45 


Corre- 
sponding 
Coefficient 
of  Variation 

0 

0.033 
O.O67 
0.100 

0.136 

0.166 
0.197 
0.230 
0.262 
0.292 

0.324 

0.351 
0.381 
0.409 
0.436 

0.462 
0.490 
0.517 

0.  544 

0.570 

0.596 
0.620 
0.643 
0.667 
0.691 

0.713 
0.734 

0.755 
0.776 
0.796 

0.818 
0.857 
0.895 
0.930 
0.966 

1.000 
1.081 
1.155 
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To  use  the  log-probability  method  in  predicting  flood  frequencies,  mean, 

standard  deviation  and  coefficient  of  skew  are  computed.    Values  of  K  from 

Table  A-II  are  then  used  to  determine  peak  discharges  corresponding  to  the 

desired  probabilities.    The  resulting  points  are  plotted  on  log-probability 

paper  to  define  the  discharge -re turn  frequency  curve.    Chow  suggests  that  if 

the  predicted  line  does  not  satisfactorily  fix  the  recorded  data,  new  values 

of  C    can  be  arbitrarily  assumed,  and  K  values  re-determined.     Lower  values 
s 

of  C    result  in  curves  that  are  increasingly  concave  upward,  while  higher 
s 

values  of  C    produce  curves  that  are  increasingly  concave  downard. 
v 

Log-Pearson  Type  III  Distribution;     The  requirement  that  all  U.S.  federal 
agencies  are  to  use  this  method  for  predicting  flood  frequencies  (Hydrology 
Committee,  Water  Resources  Council,  1967)  is  designed  to  provide  a  common  basis 
for  evaluating  proposed  water  resource  projects;  its  use  does  not  necessarily 
imply  that  it  is  the  correct  distribution  which  applies  universally  to  all 
hydrologic  data. 

The  basis  for  the  method  was  advanced  by  Pearson  in  1922  and  proposed 
for  use  in  hydrologic  studies  by  Foster  (192^) .    Actually,  Pearson 
proposed  12  different  distributions,  but  only  his  Type  III  has  received 
seru'orrss  consideration.    Beard,  (1962)  described  the  logarithmic  trans- 
formations which  produce  the  distribution  now  termed  Log-Pearson  Type  III. 

In  using  this  distribution  the  mean,  standard  deviation  and  coefficient 
of  skew  are  computed  as  before,  but  using  logarithms  of  recorded  discharges 
instead  of  the  recorded  values  themselves.      For  the  computed  coefficient  of 
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skew  and  desired  probabilities,  values  of  the  frequency  factor  K  are  deter- 
mined from  Table  A-III,  and  these  values  used  to  determine  log  Q  in 


Log  Q  =  M  +  kS  (A-lh) 

where  M  is  mean  of  logarithms  of  recorded  discharges 

S  is  standard  deviation  of  logarithms  of  recorded  discharges . 

The  probability  values  in  Table  A-III  are  only  valid  for  infinitely  long 
periods  of  record.    It  is  necessary  therefore  to  adjust  the  probabilities 
depending  on  the  length  of  record,  using  Table  A-TV. 

Extreme  Value  Distribution:     Fisher  and  Tippett  (1928)  showed 
that  if  one  selected  the  largest  event  from  each  of  many  large  samples,  the 
distribution  of  these  extreme  events  was  independent  of  the  original  distri- 
bution and  conformed  to  a  limiting  function.    In  19^1  Gumbel  suggested  fitting 
the  extreme  value  distribution  to  flood  data,  because  the  annual  flood  could 
be  assumed  to  be  the  largest  of  a  sample  of  365  possible  values  each  year. 
The  resulting  techniques  are  referred  to  as  the  "extreme  value"  or  "Gumbel" 
distribution.' 

Fisher-Tippett  argued  that  the  limiting  distribution  was  of  the  form  of 
a  double-exponential, 

-y 

P    =  e~e  (A-15) 
x 

where  P    =  probability  that  a  given  value  of  x  (peak  discharge,  for  instance) 
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TABLE  A-  IV:     TABLE  OF  PR  VERSUS  P^  IN  PERCENT 
For  use  with  samples  drawn  from  a  normal  population 
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having  small  skew  coefficients. 

Pearson  Type 

III  distributions 
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will  not  be  equalled  or  exceeded. 

e  =  base  of  naperian  logarithms 

y  =  a  function  of  x  (referred  to  as  a  reduced  variate) 
Depending  on  the  relationship  that  exists  between  x  and  y ,  the  distribution 
was  categorized  by  Fishe r-Tippett  into  three  types. 

Type  I  -  for  unbounded  data  ranging  from  -  *9  to  +  po 

Type  II  -  data  with  a  lower  bound 

Type  III  -  data  with  an  upper  bound 

For  the  asymptotic  distribution  of  Type  I,  Gumbel  (196k)  developed 
the  equation  relating  the  variable  x  to  its  reduced  variate  y ,  as 


(x  -  u) . 


(A-16) 


In  this  equation  u  is  a  "location  parameter"  given  by 


y. 


n 


u  =  x 


(A-17) 


and  oC  is  a  scale  parameter  or  measure  of  dispersion  given  by 


(A-18) 


Here    X    =  mean  of  the  sample 


-  expected  mean  of  the  population 


Sx  =  standard  deviation  of  the  sample 


(?  n  =  expected  standard  deviation  of  the  population 
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Eliminating  u  and  oC  in  equation  A-l6  gives 


(A-19) 


Gumbel  developed  formulas  for  obtaining  y    and  CP    (functions  of  the  length 

r  n  n 

of  record),  and  provided  tables  with  values  of  y^  and  (p^  for  N  (number  of 
extreme  events)  ranging  from  15  to  1000.     (See  Table  A-V) . 


In  using  the  "Gumbel"  distribution  for  determining  flood  frequencies  then, 
the  steps  are  as  follows: 

1.  Compute  sample  mean  and  standard  deviation. 

2.  Depending  on  number  of  years  of  record,  use  Table  A-V  to  obtain 
expected  mean  and  expected  standard  deviation. 

3.  For  desired  recurrence  intervals,  use  Table  A-VT  to  obtain  correspond- 
ing values  of  y. 

k.      For  each  value  of  y  solve  equation  A-19  (repeated  below  in  terms  of 
discharge)  to  obtain  required  peak  discharge  corresponding  to  each  recurrence 


A  computer  solution  of  equation  A-15,  P  55  1  -e         is  shown  in  Table 
A-VI,  where  y  values  are  tabulated  against  recurrence  interval  (t    =  ^-) . 


interval. 


3 


Q  =  Q+^(y-yn) 


(A-20) 
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TABLE  A-V:     EXPECTED  MEAN  y  ,  AND  EXPECTED  STANDARD  DEVIATION^. 
AS  FUNCTIONS  OF  THE  NUMBER  OF  EXTREMES  N 


Number 

Expected 

Expected 

of 

Mean 

Standard 

Extremes  j 

Deviation 

N 

y 

JN 

^N 

10 

0.50213 

0.97831 

15 

0.51284 

1.02057 

20 

0.52355 

1.06283 

25 

0.53086 

1.091^5 

30 

0.53622 

1.11238 

35 

0. 54034 

1. 12847 

40 

0. 54362 

l. 14132 

45 

0. 54630 

1.15185 

50 

0.5^85^ 

1. 16066 

6o 

0.55208 

1.17467 

70 

0.55^77 

1.18536 

8o 

0. 55688 

1.19382 

90 

O.5586O 

1.20073 

100 

0.56002 

1.20649 

150 

0. 5646i 

1.22534 

200 

0.56715 

1.23598 

250 

0.56878 

1.24292 

300 

0.56993 

1.24786 

J 1  nn 
4-UU 

U .  ^  f  14-4- 

1 .  djuf^)\j 

500 

0.57240 

1.25880 

750 

0.57377 

1.26506 

1000 

0.57^50 

1.26851 
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TABLE  A- VI:     PLOTTING  COORDINATES  FOR  PROBABILITY  PAPER 


y 


y 


1.01 

-1.5293 

1.02 

-1.3691 

1.03 

-l. 2630 

1  oh 

-1  1811 

1  OS 

—  j- .  iijj 

1. 10 

-0.87^5 

1.15 

-0.7114 

1.  20 

-0. 5031 

1  ?S 

-0  47^8 

1  ^0 

-0  ^827 

w .  jut:  [ 

1.40 

-0.2253 

1. 50 

-0.0940 

1.  60 

0.0193 

1  70 

0  1 1  QS 

1.80 

0  2095 

1.90 

0.2914 

2.00 

0.3665 

2 . 20 

0. 5006 

2.40 

0  61 80 

?  60 

2.80 

0.8168 

3-00 

0.9027 

3.  20 

0. 9olo 

3.40 

1.0546 

3.60 

1. 1226 

3.80 

1.1861 

4.00 

1.2458 

4 .  50 

1.3810 

5 . 00 

5  50 

1  6060 

6.00 

1.7019 

6.50 

1.7894 

7.00 

1.8698 

7.50 

1.9442 

8.00 

2.0134 

9.00 

2.1389 

10.00 

2.2503 

11.00 

2.3506 

12.00 

2.4417 

±j  .  UU 

2. 5251 

i4.oo 

2.6022 

15.00 

2.6737 

20.  00 

2.9701 

25.00 

3.1985 

^n  nn 

35.00 

3. 54o8 

40.00 

3.6762 

45.00 

3.7954 

50.00 

3.9019 

f^n  nn 
ou .  UU 

70.00 

4.2413 

80.00 

4.3757 

90.00 

4.4942 

100 . 00 

4.6001 

1  on  on 

-LciU  .  UU 

i4o.oo 

4.9380 

160.00 

5.0720 

180. 00 

5. 1901 

200.00 

5.2958 

odd  nn 

C-C\J  .  UU 

P  •  J?J-J 

240.00 

5.4785 

260.00 

5- 5587 

280.00 

5.6330 

300.00 

5.7021 

qcn  nn 
opu . UU 

P . OpDp 

4oo.oo 

5 • 9902 

450.00 

6. 1081 

500.00 

6.2136 

600.00 

6.3961 

700.00 

6.5504 

800.00 

6.6839 

900.00 

6.8018 

1000.00 

6.9072 
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Chow  has  shown  that  the  extreme  value  distribution  is  really  essentially 
a  special  case  of  the  log-probability  law  where  the  skew  coefficient  is  1.139, 
and  coefficient  of  variation  is  O.36U. 

The  "Gumbel"  distribution  has  been  widely  publicized  in  hydrology  text- 
books and  is  frequently  used  today. 

In  order  to  minimize  bias  involved  in  the  prediction  of  extremes, 
Lieblein  (195*0  suggested  a  method  to  fix  the  extreme  value  distribution 
which  preserves  the  original  time  order  of  occurrence  of  events. 

I.  I.  Gringorten  (1962)  suggested  that  no  single  distribution  can  be 
expected  to  fit  all  meteorological  data  universally.    He  developed  seven 
"type  curves", to  which  recorded  data  are  compared,  and  the  distribution 
best  fitting  the  data  is  thereby  selected.     Four  of  the  curves  are  essentially 
"normal  distribution"  curves.    Curve  1  is  the  basic  normal  distribution; 
curve  2  is  the  distribution  of  the  largest  of  10  members  taken  independently 
from  a  normal  distribution;  3  and  k  are  the  distributions  of  the  largest 
of  100  and  1000  members  respectively,  taken  independently  from  a  normal 
distribution.     Curve  5  is  the  double  exponential  distribution,  whiie  6  and 
T  are  versions  of  the  Fisher-Tippett  Type  II  distribution. 

The  fundamental  equations  governing  the  family  of  Gringorten  curves  can 
be  stated  as 

x  -  x  +  (T   y  (A-21 ) 

x  r 
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where 


yr  -  (y  -  y)A'y 


(A-22) 


In  these  equations  x 


estimate  of  the  variable 


x 


true  mean  of  x 


true  standard  deviation  of  x 


y 


standardized  value  of  idealized  variate 


y 


true  mean  of  y 


true  standard  deviation  of  y. 


Figure  A-h  shows  the  type  curves,  while  Table  A -VI I  gives  descriptions  of  the 
curves  and  their  characteristic  features. 

Hydrologic  Simulation 

In  the  last  few  years  the  rapid  development  of  mathematical  statistics 
and  the  digital  computer  have  made  it  possible  to  generate  hypothetical  sequences 
of  events  which  have  the  same  probability  characteristics  as  the  historic 
record,  and  thereby  permit  the  consideration  of  many  possible  combinations  of 
output.     Probably  the  best-known  simulation  analysis  to  date  is  that  the 
Harvard  Program,  as  reported  by  Maass,  et  al  (1966).    With  32  years  of 
record  available  for  the  Clearwater  River  in  Idaho,  a  500-year  streamflow 
sequence  was  generated  by  mathematical  synthesis. 

Several  stochastic  models  have  been  proposed  for  simulation,  falling 
into  one  of  two  general  categories.     If  the  individual  items  of  data  are 
considered  statistically  independent  (that  is,  the  state  of  the  system  at  a 
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given  time     is    not  dependent  upon  previous  states)  the  process  is  said  to 
be  "non-Markovian , "  and  the  usual  treatment  is  application  of  Monte  Carlo 
methods.     If  there  are  carry-over  effects  (the  streamflow  at  one  time,  for 
instance,  being  dependent  upon  the  streamflow  at  an  earlier  time)  the  Monte 
Carlo  methods  cannot  be  applied  directly,  and  use  of  transition  probability 
theory  is  usually  used  instead.     Transition  probability  is  developed  under 
the  general  theory  of  Markov  processes . 

Monte  Carlo  methods  -  After  studying  the  characteristics  of  the  distribution 
of  recorded  data  an  appropriate  sampling  technique,  involving  the  use  of  tables 
of  random  numbers,  is  used  to  produce  artificial  sequences  having  the  same 
stochastic  structure  as  the  record.    Fie  ring  (196?)  used  Monte  Carlo 
techniques  in  studying  the  Lehigh  River  basin  in  Pennsylvania. 


Markov  Processes  -  A  sample  expression  of  a  Markov  model  is 

yt  =  p^-i  +vf  s  (A_23) 

where  y    -  runoff  at  time  t 

y  f  ±  ~  runo^f  a-t  "time  T-  1 
p^  =  constant 

r\    =  constant  dependent  on  the  distribution  of  the  recorded  runoff  dat,a 
E  =  a  random  variable 

The  coefficients  must  be  properly  chosen  so  that  the  simulated  sequences 
will  have  the  same  stochastic  properties  as  the  recorded  series .  Additional 
sophistication  may  be  obtained  by  adding  terms  (i.e.,  p    yr_P»  PJ-_m  etc.). 
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Amorocho  and  Hart  suggest  that  there  is  no  clearcut  way  to  decide, 
in  general,  whether  to  assume  no  time  dependence  and  use  Monte  Carlo  methods 
directly,  or  whether  to  assume  simple  or  complex  time  dependence. 

SUMMARY 

The  purpose  of  this  chapter  has  been  to  describe  the  current  techniques 
available  for  determining  peak  discharge  and  associated  frequency.  The 
parametric  methods  permit  calculation  of  discharge  peak,  but  do  not  suggest 
a  return  frequency  or  probability.     In  using  parametric  methods  a  frequency 
may  be  assigned  to  a  causative  storm,  and  this  will  permit  indirectly  the 
determination  of  peak  flow  frequency.     Stochastic  methods  use  the  statistical 
characteristics  of  the  recorded  data,  precipitation,  or  streamflow. 
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APPENDIX  B 


1.  Map  of  Montana  Showing  Instrumented  Watersheds 

2.  Maps  of  Instrumented  Watersheds 

a.  Bacon  Creek 

b.  Duck  Creek 

c.  Hump  Creek 

d.  Lone  Man  Coulee 

3.  Aerial  Photographs  of  Instrumented  Watersheds 

a.  Bacon  Creek  looking  Northwest 

b.  Bacon  Creek  looking  Southeast 

c .  Duck  Creek  looking  South 

d.  Duck  Creek  looking  North 

e.  Hump  Creek  looking  South 

f.  Hump  Creek  looking  North 

g.  Lone  Man  Coulee  looking  West 

h.  Lone  Man  Coulee  looking  East 
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FIGURE  B-l:     Bacon  Creek,  looking  Northwest 


FIGURE  B-2:     Bacon  Creek,  looking  Southeast 


FIGURE  B-k:     Duck  Creek,  looking  North 


FIGURE  B-5:     Hump  Creek,  looking  South 


FIGURE  B-6:     Hump  Creek,  looking  North 


FIGURE  B-7:     Lone  Man  Coulee,  looking  West 


FIGURE  B-8:     Lone  Man  Coulee,  looking  East 
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Watershed  Characteristics 

a.  Area 

b.  Shape 

c.  Compactness  Coefficient 

d.  Orientation 

e.  Elevation 

f.  Stream  Channel  Length 

g.  Stream  Channel  Slope 

h.  Overland  Slope 

i .  Number  of  Streams 
j .  Drainage  Density 
k.  Stream  Density 

1.  Land  Use 

m.  Infiltration  Index 

n.  Ponds 
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APPENDIX  D 


SOIL  CHARACTERISTICS 

a.  Grain  Size  Analysis 

b.  Atterberg  Limits 

c.  Soil  Survey  Data 

d.  Water  Intake  Rates 
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BACON  CREEK  WATERSHED 


SOIL  SURVEY  DATA 


Publication 
Symbol 

Aa 

Ab 
Ac 

Ad 

Ba 
Bb 

Be 
Bd 
Be 
Bf 
Ca 
Cb 
Cd 
Ce 
Ma 
Mb 
Mc 
Md 
Me 
Mf 
Mg 
Ta 


Hydrologic 
Soil  Group 


Soil  Name 


B 


B 


D 


Alluvial  land  clay  and  saline, 

2-1+%  slope 
Ashuelot  gravelly  loam,  0-2%  slope 
Ashuelot-Musselshell  complex, 

0-2%  slope  - 
Ashuelot-Musselshell  complex, 

2-k%  slope 
Bercail  clay  loam,  l-h%  slopes 
Bercail  clay  loam,  shale  substratum, 

k-10%  slopes 
Bercail-Thebo  clays,  k-10%  slopes 
Bercail-Thebo  clays,  10-15%  slopes 
Bercail-Thebo  clays,  15-1*5%  slopes 
Bercail  clay  loam,  saline  1+-1Q%  slopes 
Crago  loam,  1-1+%  slope 
Crago  loam,  l+-8%  slope 
Crago  cobbly  loam,  8-15%  slope 
Crago-Bercail  complex,  15-1+5%  slope 
Musselshell  gravelly  loam,  0-2%  slope 
Musselshell  gravelly  loam,  2-k%  slope 
Musselshell  gravelly  loam,  l+-8%  slope 
Musselshell  gravelly  loam,  8-15%  slope 
Musselshell-Crago  complex,  0-2%  slope 
Musselshell-Crago  complex,  2-8%  slope 
Musselshell-Crago  complex,  8-15%  slope 
Thebo  clay,  1+-2Q%  slopes 

Total  Acres 


Approximate 
Acres 


280 
258 

1169 

37^ 
1207 

380 
1000 

19^ 
11+79 
99 

216 

89 
kkl 

3099 
2366 
989 
872 
116 
2616 

75^ 
178 
1+88 
18661+ 


1.5 
1.1+ 

6.3 


-  D-2  - 


1 


] 

DUCK  CREEK 
SOIL  SURVEY  DATA 


] 


.ication 

Hydrologic 

Approximate 

=^ibol 

Soil  Group 

Soil  Name 

Acres 

% 

j 

"a 

B 

Alluvial  land,  loamy  0  to  h%  slopes 

1337 

3.85 

Ba 
-i 

B 

Briley  loam,  2  to  h%  slopes 

9 

0.03 

}b 

B 

Briley  gravelly  loam,  2  to  h%  slopes 

33 

0.10 

B 

Brockway  silty  clay  loam,  0  to  2%  slopes 

275 

0.79 

Bd 

B 

Brockway  silty  clay  loam,  2  to  k%  slopes 

h5 

0.13 

Ta 

B 

Cherry  silty  clay  loam,  0  to  2%  slopes 

12k0 

3.57 

> 

B 

Cherry  silty  clay  loam,  2  to  h%  slopes 

87 

0.25 

Cc 

B 

Cherry  silt  loam,  0  to  2%  slopes 

^5 

0.13 

"T  d 

B 

Cherry  silty  loam,  2  to  k%  slopes 

267 

0.77 

3e 

B 

Chinook  fine  sandy  loam,  1  to  6%  slopes 

11*6 

0.^2 

"fcg 

C 

Cabba  complex,  15  to  k0%  slopes 

2720 

7.83 

Ch 

B/C 

Cabba- Brandenburg  complex,  8  to  ^5%  slopes 

70 

0.20 

B/C 

Cabba-McGuire  complex,  k  to  8%  slopes 

165 

0.U8 

J3J 

B/C 

Cabba-McGuire  complex,  8  to  20%  slopes 

U05 

1.17 

Ck 

B 

Cabba-Pennel  loams ,  8  to  20%  slopes 

3721+ 

10.72 

I1 

C/D 

Cabba-shale  outcrop  complex,  15%  plus  slopes 

11.67 

w 

B 

Farnuf  loam,  0  to  2%  slopes 

813 

2.3h 

Fb 

B 

Farnuf  loam,  2  to  h%  slopes 

591 

1.70 

.Ha 

B 

Harkan  gravelly  loam,  6  to  35%  slopes 

hlk 

1.19 

ffb 

B/C 

Harkan-Cabba  complex,  15  to  h0%  slopes 

880 

2.53 

•He 

B/C 

Harkan  and  Cabba  soils,  6  to  15%  slopes 

1380 

3.97 

La 

C 

Lambert  and  Lothair  soils,  15%  and  steeper  slopes 

590 

1.70 

Ma 

A 

McGuire  fine  sandy  loam,  h  to  8%  slopes 

101 

0.29 

A 

I'ituuii  c  line  sdnciy  ±oajn,  o  uo  d.\j/o  S-Lopes 

4  O  -J 

Mc 

A 

McGuire  loamy  fine  sand,  k  to  20%  slopes 

322 

0.93 

■Md 

B 

Macay  loam,  0  to  2%  slopes 

1659 

^.78 

Me 

B 

Macay  loam,  2  to  h%  slopes 

318 

0.92 

Pa 

B 

Pennel-Macay  loams ,  2  to  k%  slopes 

96 

0.28 

Pb 

B 

Pennel-Macay  loams,  k  to  8%  slopes 

1621+ 

U.  68 

Pc 

B 

Pennel-Cabba  loams,  k  to  8%  slopes 

715^ 

20.60 

'Ra 

B 

Remuda  fine  sandy  loam,  2  to  8%  slopes 

698 

2.01 

Rb 

B 

Remuda-McGuire  complex ,  k  to  8%  slopes 

1620 

h.6l 

Rc 

B 

Romine  fine  sandy  loam,  0  to  2%  slopes 

70 

0.20 

Ya 

B/C 

Yawdim-Cabba  soils  k  to  35%  slopes 

199 

0.57 

Total 

3^,721* 

100.00 

-  D-3  - 


HUMP  CREEK  WATERSHED 


SOIL  SURVEY  DATA 

Publication  Hydrologic  Approximate 

Symbol  Soil  Group   Soil  Name   Acres  %  j 

Aa  C  Absarokee-Amherst ,  2-8  per  cent  slope  259  5°  8 

Ab  Absarokee-Amherst,  8-20  per  cent  slope  132  2.9 

Ac  Absarokee-Amherst,  stony  clay  loam  k-Q%  slope  3^8  7.8 

Ad  C  Amherst-Castner ,  stony  loams  8-15%  slope  720  16.1 

Ba  C  Basket  clay  loam  complex,  2-8%  slope  118  2.6 

Bb  C/D  Beaudry-Absarokee  complex  U-8%  slope  60  1.3f 

Ae  C  Absarokee  clay  loam,  clay  substratum  4-8%  slope  19  0.h\ 

Ma  C/D  Midway-Castner  association,  k-10%  slope  6k0  lk.3 

Mb  C/D  Midway-Castner  association,  10-30%  slope  509  11. 3i 

Mc  D  Midway-Rock  outcrop  complex,  15-50%  slope  1565  35.0 

Wa  C  Work  clay  loam,  4-8%  slope  71  1.6 

Wb  Work  clay  loam,  Q-l6%  slope  _3£  0.9, 

Total  acres  kkQO  100.  of 


LONEMAN  COULEE  WATERSHED 


SOIL  SURVEY  DATA 


folic at ion        Hydrologic  Approximate 
Symbol  Soil  Group   Soil  Name   Acres  % 


Aa 

Alluvial  land,  poorly  drained  and  saline, 

u  XiO  Lt/o  S-Lopes 

DOU 

1  •  -> 

Ca 

C 

Cabba  loam  and  clay  loam,  8-20%  slopes 

281 

3.0 

Cb 

Cabba  loam-Rockland  complex,  15-^0$  slopes 

107 

1.1 

Ka 

B 

Kiev  loam,  0  to  2%  slopes 

2283 

2h.k 

Kb 

Kiev  loam,  2  to  k%  slopes 

IU98 

16.0 
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APPENDIX  E 


HYDROLOGIC  DATA 

a.  Hydrographs 

b.  Antecedent  Meteorological  Conditions 

c.  Monthly  Precipitation  Totals 

d.  Snow  Survey  Measurements 
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APPENDIX  F 


VARIABLES  FOR  MULTIVARIATE  STATISTICAL  STUDIES 

a.  Description 

b .  Value  s 


APPENDIX  F 


Descriptions  of  Watershed  Parameters  and  Storm  Variables  Used  in 

Mult i-Vari ate  Statistical  Studies 

WATERSHED  PARAMETERS 


1.  A      :     Area  of  watershed  (sq  mi). 

2.  SHP  :     Shape  of  watershed  ( dimensionless )  from  the  ratio  of  the 

length  of  the  longest  line  that  could  be  passed  through  the 
center-of-area  of  the  watershed  to  the  length  of  a  perpend- 
icular through  this  point.     Center-of-area  found  from  the 
balance  point  of  a  cut-out  of  the  watershed  from  SCS  soil 


3.    AZ     :     Azimuth  of  main  channel  (deg)  measured  clock-wise  from  true 
North  when  referred  to  the  water -stage  recorder.  Mean 
azimuth  of  straightline  sections  of  main  channel  computed  from 

(AF.  L.  ) 
^— -       i  i 

AZ  - 


Z 


where 


AZ^  =  azimuth  of  each  straight-line  segment,  starting  at 
gaging  station  and  ending  at  closest  point  on  basin 
divide  to  end  of  main  channel  on  SCS  soil  map 


-  F-l  - 


L. 
i 


=  length  of  each  straight  line  segment. 


b.    ELEV:     Elevation  of  watershed  (ft)  computed  from  an  average  of  the 

elevations  at  0.2  and  0.8  of  the  main  channel  straight-line-segment 
length  from  the  gaging  station  to  the  nearest  point  of  basin 
divide  from  head  of  main  channel. 

5*     °NDS:     M^imum  ground  slope,  (ft/ft)  computed  by  dividing  the  distance  from 
the  head  of  the  blue,  main-channel  line  on  USGS  Quadrangle  map  to 
the  nearest  point  on  the  basin  divide,  by  the  change  in  elevation 
between  these  points. 

6*     GNDL:     Overland  distance  (mi)  of  flow  obtained  from  distance  from  head 

of  blue,  main-channel  line  on  USGS  Quadrangle  map  to  nearest  point 
on  basin  divide. 

T"    FREQ:     ^ream  frequency,  1/sq  mi),  or  number  of  stream  segments  per  unit 

area.     Obtained  by  dividing  the  total  number  of  stream  segments  of 
all  orders  by  the  area  of  the  watershed. 

8*    S      :     MalE  channel  slope_  (ft/ft)  computed  as  the  slope  of  a  straight 
line  drawn  through  the  main  channel  straight-line  -segment 
profile.    The  straight  line  is  fitted  to  the  profile  in  such 
a  way  as  to  give  the  same  area  under  the  line  as  the  total  area 
under  the  profile,  using  equal  base  lengths. 


-  F-2  - 


Main  channel  meander  length  (mi),  including  distance  from  end 


10.  USE  : 


11.  INFR : 


12 .    POND : 


of  blue  line  to  nearest  point  on  basin  divide. 

Land  use  ratio  ( dimensionless )  obtained  from  a  ratio  of  the 
vegetated  (including  forests  for  Hump  Creek)  area  to  the  bare- 
cultivated  (summer  fallow)  area.    Relative  areas  were  deter- 
mined subjectively  from  each  section  in  the  aerial  photographs . 

Infiltration  rate  (in./hr)  of  water  shed,  weighted  for  rates 
different  soil  types.    Computed  from 


I.  =  average  measured  infiltration  rate 
i 

on  each  soil  type  of  watershed.     If  not 
measured,  soil  was  classified  by  SCS  method  (type 
A,  B,  C,  or  D),  and  the  average  of  all  watersheds  was 
applied  to  any  soil  types  which  were  not  tested 

=  area  of  each  soil  type  on  watershed  from  planimeter  of 
SCS  soil  map. 

Distance-weighted  per  cent  of  total  area  in  ponds  and  reservoirs 
{%)  computed  from 


-  F-3  - 


( L.  AP.) 


i 


POND  = 


A2 


(100) 


where 


L.  =  straight-line-segment  distance  along  main  channel 
from  gaging  station  to  each  pond  on  SCS  soil  map 

AP.  =  full  surface  area  of  each  pond  or  reservoir 

i 

A      =  area  of  watershed. 


STORM  VARIABLES 


13.  I      :     Intensity  of  precipitation  (in./hr)  computed  as  mean  of  all 


hourly  intensities  from  all  recording  precipitation  stations  from 
beginning  of  storm  on  the  first  recording  station  to  hour  of  last 
runoff  from  hydrograph . 


ik .  ISD  :     Standard  deviation  of  precipitation  intensties  used  in  variable 


13,  computed  from 


1 


n-I 


n 


where 


n  =  total  number  of  hourly  intensities  recorded  for  all 
recording  precipitation  stations. 


15-  D       :     Duration  of  storm  (hr)  computed  as  an  average  time  to  the 

center  of  area  of  rainfall  hyetographs  for  all  recording  prec- 
ipitation stations s  computed  from, 


5"  (  DCA.  A.  ) 


D  = 


Where 


(  I.   T.  ) 

1  I 


i. 


DCA,  -  time  in  hours  from  beginning  of  each  rainfall 
x 

hyetograph  to  the  center  of  area  of  each 
hyetograph,  where  hyetograph  ends  with  last 
precipitation  before  hour  of  end  of  runoff 


A    =      Thiessen  area  of  each  recording  precipitation 
station , 


T.  =      time  in  hours  from  beginning  of  hyetograph  to 

i. 

I , ,  where 


I .  =      hourly  precipitation  at  the  respective  recording 
s  tat ion • 

16 .  TDF  :     Time  distribution  factor  (hr)  of  precipitation  over  watershed, 
computed  as , 


-  F-5  - 


TDF  -  Dl  -  D 


where 


and 


y  (dca  a.) 

Dl  =   ^  

A1 


DCA,  =  time  in  hours  from  beginning  of  first  appearing 

J 

rainfall  hyetograph  to  the  center  of  area  of  each 
hye to graph 


A    =  Thiessen  area  of  each  recording  precipitation 
J 

station . 


IT.  TPCP:     Total  precipitation  (in.)  computed  as  an  average  Thiessen- 
weighted  precipitation  from  all  precipitation  stations, 


y  (a.  pj 

TPCP  = 


where 


P^  =  total  storm  precipitation  at  each  recording  and  non- 
recording  precipitation  station  from  beginning 
of  storm  to  end  of  runoff 
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A.  =  Thiessen  area  of  each  station  from  planimeter  of 
1 

SCS  soil  map. 


18.    API  :     Fourteen-day  antecedent  precipitation  index  (in.)  computed 
from 


where 


API.  =  k  (API.  J  +  (P. ) 
1  l-l  1 


API^  =  antecedent  precipitation  index  for  the  ith  day 

after  the  fourteenth  day  before  the  first  day  of 
the  storm 


k  =  0.78,  a  reduction  constant  to  indicate  the 

evaporation-transpiration  losses  from  day  to  day. 
(Usually  ranges  from  0.80  to  0.95  (Chow,  p.  25-102), 
and  small  value  used  here  for  relatively  dry  water- 
sheds . ) 


P.  =  Thiessen -weigh ted  average  daily  precipitation  for  each 

X 

station 


API11+  =  API  for  first  day  of  storm. 
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19.  SOLM:     Average  water  content  {%)  of  soil  for  all  stations  for  3,  9,  and 

1 8-inch  depths,  at  beginning  of  first  hyetograph. 

20.  WDIR:    Predominant  wind  direction  during  storm  period  (dimensionless ) , 

obtained  from  weather  station  data,  and  given  the  numbers  1, 
2,  3,   . .  .  ,  8  to  represent  a  wind  from  the  N,  NW,  W,  SW,  S,  SE, 
E,  NE,  respectively. 

21.  WEEK:     Week  of  the  year  (dimensionless)  in  which  the  date  of  the 

peak  discharge  rate  fell.  Weeks  were  numbered  from  the  first  week 
in  Jan, 

22.  AIRT:     Average  Thiessen-wei ghted  air  temperature  (deg  F),  h  ft  above 

ground,  during  storm  period.     Given  by 

£(T.  A.) 
AIRT  =  C 


where 


T\  =  mean  hourly  temperature  at  a  station  during  the  time 

from  beginning  of  first  precipitation  to  end  of  runoff 
hydrograph 


A.  =  Theissen  area  of  each  weather  station  on  the  watershed 

i 
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23.  ATSD:     Average  Thi esse n-weigh ted  standard  deviation  of  air  temperatures 
(deg  F),  during  the  storm  period.     Given  by 


V2 

2  1      !  <Tm  2  -  v^  i 


(£  T 

i  n-1      ^  i 


n 


ATSD  = 


91  ST.  A.) 
^  ii 


Z  Ai 


where 


ST^  =  standard  deviation  of  hourly  air  temperatures  at 
a  weather  station  during  the  storm  period 


n  -  number  of  hourly  air  temperatures  considered 


T.  =  hourly  air  temperature  at  k  ft  above  ground  at  a 
i 

station 


A.  =  Thiessen  area  of  each  weather  station  on  the 
i 

watershed. 


2k,  WVEL:     Mean  wind  velocity  (mph)  during  the  storm  period. 

25.  WVSD:     Standard  deviation  (mph)  of  hourly  wind  velocities  during  the 
storm  period. 


-  F-9  - 


26.  SOLT :     Average  Thi ess en-weighted  soil  temperature  (degF),  3  in.  below 

surface,  during  storm  period.     Computed  exactly  as  air 
temperature  mean,  only  using  four  readings  per  day  at 
6-hr  intervals,  beginning  at  3:00  a.m. 

27.  STSD:     Average  Thiessen-weighted  standard  deviation  of  soil  temperatures 

(deg  F),  3  in.  below  surface,  during  storm  period.     Computed  in 
the  same  manner  as  the  air  temperature  standard  deviation. 

28.  DEGD:    Average  Thiessen-weighted  Degree-Days  (deg  F)  for  a  1^-day 

period  prior  to  the  beginning  of  runoff.  Degree-Days  for  each  day 
was  maximum  daily  temperature,  and  l^+-day  average  was 

ST  (DD.  A.  ) 

DEGD  =  - 


where 


DD.  -  total  Degree-Days  for  a  weather  station  for  a 
l 

1^4-day  period  before  the  first  day  of  runoff 


A.  =  Thiessen  area  of  each  weather  station  on  the 
i 

watershed. 


29.  SWEQ:     Snow-water  equivalent  (in.),  or  the  volume  of  water  in  the  snow 
pack  on  the  watershed  on  the  date  of  the  beginning  of  runoff. 
Computed  from  Thiessen-weighted  volume  from  each  snow  course. 


-  F-10  - 


Loss  or  gain  in  volume  from  date  of  snow  survey  to  date  of  runoff 
was  computed  as  the  difference  between  accumulated  precipitation 
and  accumulated  runoff.    This  was  added  to  the  Thiessen-weighted 
volume  on  the  watershed  at  the  time  of  the  last  snow  survey, 
giving  the  volume  on  the  date  of  beginning  of  runoff. 
Thiessen  areas  of  snow  courses  were  used  to  give  the  volume  on 
the  data  of  the  survey, 


V 


where 


VC.  =  snow-water  equivalent  at  a  snow  course 
1  ^ 


A. 


Thiessen  area  of  each  snow  course  on  the  watershed. 


i 
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APPENDIX  G 


COMPARISON  OF  STOCHASTIC  PEAK  FLOW  METHODS 


STATISTICAL  ANALYSIS  OF  SEVEN  STOCHASTIC  METHODS 

Bacon  Creek  Watershed 

Duck  Creek  Watershed 

East  Fork  Duck  Creek  Watershed 

Hump  Creek  Watershed 

Lone  Man  Coulee  Watershed 
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APPENDIX  J 


EFFECT  OF  RECORD  LENGTH  ON  PEAK  FLOW  PREDICTIONS 

Grasshopper  Creek 
Pryor  Creek 
Burnt  Fork  Creek 
McDonald  Creek 


TABLE  J-l  EFFECT  OF  RECORD  LENGTH  ON  PEAK  FLOW  PREDICTIONS  AT  GRASSHOPPER 
CREEK     (348  square  miles,  33  years  record) 
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1st  16  years  record 
2nd  17  years  record 
total  33  year  record 


Mean  Annual  Peak,  cfs 

59k-7 
^75.2 


Standard  Deviation,  cfs 
192.9 
397.7 
33k.9 


Gumbel  Method: 

1st  l6  years  record 
total  33  year  record 

Log-r Normal  Method: 

1st  l6  years  record 
total  33  year  record 

Log- Pearson  Method: 

1st  l6  years  record 
total  33  year  record 


I.Q5  yr 

60.9  cfs 
126.9 

1 . 05  yr 
107.1  cfs 
126.3  cfs 


6i.5  yr 

234.1  cfs 

286.3 

Q 

^2.0  yr 


330.9  cfs 
391.5 


0, 


J2.0  yr 
296.8  cfs 
364.6 


^200  yr 
1244.5  cfs 
1895.2 

J  Q 

^5  yr  ^100  yr 

502.6  cfs  869.1  cfs 

669.4  1690.9 

^2.5  yr      Ql6  yr 
408.2  cfs  817.1  cfs 
510.2  1052.0 


-  J-l  - 


■   __p  -   Qer.  t  Order  t  

!912  930  6  3  83  -2 

!913  1130  k  5  75  13  3.46 

L9^  595  14  1*64  10  4.50 

J*1!  6^  10  2.3o  26  1.73 
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3             7'67  3  15.*00 

1st  22  years  record  ^L^^PeaK  cfs  Standard  Devlin, ,  cfs 

2nd  22  years  record  fi^HT^   

total  44  year  record  7ol  n  793-l 

f93-°  608.5 
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TABLE  J-  2(  Continued) 

EFFECT  OF  RECORD  LENGTH  ON  PEAK  FLOW  PREDICTIONS  AT  PRYOR 
CREEK     (^38  square  miles,  kk  years  record) 


Gumbel  Method: 


G1.5  yr  ^200  yr 


1st  22  years  record  508.9  cfs  2175.2  cfs 

total  kk  year  record  k^k.k  cfs  3306.9  cfs 

Log-Normal  Method:  Q  Q  ft  ft 

—    ^--05  yr  2 . 0  yr  5  100  yr 

1st  22  years  record        23^-3  cfs  657.6  cfs  956 . 8  cfs     1665. 1  cfs 

total  Ml  year  record      221.0  622.6  1097-3  3099-2 

Log- Pearson  Method:  ft,  nc  Qn  _  ft0  „  ft. „ 

— 2   ^1.05  yr  ^2  .  0  yr  3 . 3  yr        17  yr 

1st  22  years  record        223-2  cfs  668.7  cfs  9^7-0  cfs    2004  cfs 

total  kk  year  record      188.1  652.3  967-5  2262 


-  J-3  - 


TABLE  J-3  EFFECT  OF  RECORD  LENGTH  ON  PEAK  FLOW  PREDICTIONS  AT  BURNT  FORK 
CREEK     (74  square  miles,  36  years  record) 


Year  Discharge 
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_Gumbel  Method: 
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Log-Normal  Method: 

1st  18  years  record 
total  36  year  record 

Log- Pearson  Method: 
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20 

1 .  yu 

3.17 

10 

J  .  ou 

1.36 

24 

1.58 

£.  L 

-L .  OU 

9.50 

2 

18  00 

6.33 

6 

6.00 

1.12 

29 

1.29 

19.00 

1 

36.00 

1.90 

19 

2.00 

4.75 

7 

5.14 

Q 

^1.05  yr 
155.7  cfs 
170.3 

Q1.05  yr 
111.3  cfs 
146.7 


Hi -I  c 

1-5  yr 

246.2  cfs 
268.3 

^2.0  yr 
306.8  cfs 

339.3 

^2.0  yr 
278.1  cfs 
325.1 


Standard  Deviation,  cfs 

164.2   \ 

208.2 
189.3 

^200  yr 
1099.2  cfs 
1166.7 

5  yr 
448 . 0  cfs 
488.3 

3-3  yr 
371.8  cfs 
418.3 


^100  yr 

907.7  cfs 
1046.7 

17  yr 
694.9  "cfs 
720.0 


-  j-4  - 


TABLE 

J-4  EFFECT 

OF  RECORD 

LENGTH 

ON  PEAK 

FLOW  PREDICTIONS  AT 

Mcdonald 

CREEK 

(^-21  square  miles 

,  31  years  record) 

First  Half  Data 

Second 

Half  Data 

Tot  &  1 

Record  Dat; 

cf  s 

Order 

t 

Order 

t 

Order 

t 

J? 

p 

_P 

1931 

4.2 

l4 

1.07 

30 

1.07 

1932 

562 

3 

5.00 

7 

4.57 

1934 

351 

8 

1.88 

16 

2.00 

1935 

371 

7 

2.14 

15 

2.13 

1936 

220 

11 

1.36 

21 

1.52 

1937 

1+01 

6 

2.50 

13 

2.46 

1938 

328 

9 

1.67 

17 

1.88 

1939 

270 

10 

1.50 

18 

1.78 

19^0 

75 

13 

1.15 

27 

1.19 

19^1 

482 

5 

3.00 

11 

2.91 

19^2 

900 

1 

15.00 

2 

16.00 

19^3 

5^9 

4 

3.75 

8 

4.oo 

1944 

853 

2 

7.50 

4 

8.00 

1945 

180 

12 

1.25 

24 

1.33 

1953 

424 

7 

2.57 

12 

2.66 

195^ 

47 

16 

1.13 

29 

1.10 

1955 

24l 

10 

1.80 

20 

1.60 

1956 

88 

14 

1.29 

26 

1.23 

1957 

67 

15 

1.20 

28 

l.l4 

1958 

150 

13 

1.38 

25 

1.28 

1959 

400 

8 

2.25 

14 

2.29 

i960 

200 

12 

1.50 

23 

1.39 

1961 

0.1 

17 

1.06 

31 

1.03 

1962 

1+85 

6 

3.00 

10 

3.55 

1963 

250 

9 

2.00 

19 

1.68 

1964 

r44o 

l 

18.00 

1 

32.00 

1965 

800 

4 

4.50 

6 

5-33 

1966 

200 

12 

1.50 

23 

1.39 

1967 

890 

2 

9.00 

3 

10.67 

1968 

810 

3 

6.00 

5 

6.4o 

1969 

500 

5 

3.60 

9 

3.56 

Mean  Annual  Peak,  cfs 

Standard 

Deviation, 

1st  l4  years  record 

396.1 

267.8 

2nd  17  years  record 

411.3 

382.8 

total  31  year  record 

4o4.5 

327.8 

Gutnbel  Method: 

Q 

^1.5  yr 

Q 

^200  yr 

1st  l4  years  record 

233.3  cfs 

1649.7 

cfs 

total  31  year  record 

219.2  cfs 

1801.5 

cfs 

Log 

Normal  Method: 

0 

^1.05  yr 

Q 

^2.0  yr 

^5  yr 

\00  yr 

1st  l4  years  record 

=  2.9  cfs 

372.0  cfs 

610.3  cfs 

1119.2  cfs 

total  31  year  record 

-15.1 

348.8 

640.5 

1446.9 

Log 

^-Pearson  Method: 

0 

^1.05  yr 

Q 

^2.0  yr 

Q3-2  yr 

Ql6  yr 

1st  l4  years  record 

28.5  cfs 

260.4  cfs 

525.3  cfs 

2381.0  cfs 

total  31  year  record 

10.9 

215.7 

554.7 

424o.o 

